o 

(N 
o 



Study of the reaction pp ppn^ within 10 MeV above 

the threshold 

COSY-TOF collaboration 
S. AbdEl Samad% R. Bilger'*, A. Bohm^ K.-Tli. Brinkmann^ H. Clement'^ 



2 ■ S. Dsliemuchadse-^ , W. Eyricli^, D. Filges*^, H. Freiesleben'', M. Fritsch^, R. Geyer' 



5 

e 



D. Hesselbarth^, B. Jakob^, L. Karsch'', K. Kilian'^, H. Koch", J. Kress'^, E. Kuhlmann^* 
S. Marwinski^ P. Michel-^, K. MoUer^, H.R Morsch^ L. Naumann^, M. Ricllter^ 
O ■ E. Roderburg*^, M. Rogge^, A. Schamlott-^, M. Schmitz^, P. Schonmeier^, W. Schroeder*^, 

^ ; M. Schulte-Wissermann^ M. Steinke", F. Stinzing^ G.Y. Sun^ G.J. Wagner'^, 

M. Wagner^ A. Wilms'', S. Wirtli'^ 

(N 
> 



" Institut fiir Experimentalphysik, Ruhr-Universitat Bochum, D-44780 Bochum, Germany 



, Institut fiir Kern- und Teilchenphysik, Technische Universitat Dresden, D-01062 Dresden, Germany 

I ° Physikalisches Institut, Universitat Erlangen, D-91058 Erlangen, Germany 
, Physikalisches Institut, Universitat Tiibingen, D-72076 Tiibingen, Germany 

I ^ Institut fiir Kernphysik, Forschungszentrum Jiihch, D-52425 Jiihch, Germany 

' Institut fiir Kern- und Hadronenphysik, Forschungszentrum Rossendorf, D-01314 Dresden, Germany 



April 25, 2003 



Abstract 



^ ' Kinematically complete measurements of the pp ppi:^ reaction were performed for beam energies 
I in the range 292 — 298 MeV. We detected both protons in coincidence by using the large acceptance 
COSY-TOF spectrometer and an external proton beam at COSY-Jiilich; thus we measured total 
and differential cross sections and energy distributions. A strong enhancement was observed in 
the Dalitz plots due to the final state interaction between the two outgoing protons; the data are 
well reproduced by Monte Carlo simulations using a standard scattering length of ao=-7.83 fm 
and effective range of ro=2.8 fm. Our measured total cross sections are roughly 50 % larger than 
those of recent internal target experiments at lUCF and CELSIUS. This discrepancy may be due 
to the final state interaction pushing some events into the very small-angle region which is almost 
inaccessible to an internal target experiment. In the angular distributions we found only very slight 
deviations from isotropy. The analysis of our momentum distributions yielded a small fraction of 
5 lb 5% for Ps or Pp wave contributions. 

PACS numbers: 13.75.Cs, 25.10.-bs, 25.40.Ep, 29.20.Dh 
# corresponding author, e.kuhlmann@fz-juelich.de 
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I. Introduction 

The last decade saw a renaissance of near-threshold pion producing reaction studies in 
nucleon-nucleon (NN) collisions. With the advent of medium-energy accelerators in Bloom- 
ington/USA (lUCF), Uppsala/ Sweden (CELSIUS), and Jiihch/Germany (COSY) a weahh 
of high precision data on total and differential cross sections as well as polarization observ- 
ables was measured for the reactions pp ppn^ [ 1-6 ], pp pmr^ [ 7-9 ] and pp —* dir^ 
[ 10-11 ], which in turn induced a flurry of theoretical activity [ 12-18 ]. So far, good agree- 
ment is found in the theoretical description of the reactions pp dn^ and pp — > pmT~^, 
larger discrepancies are observed for the pp ppn'^ reaction. Soon after publication of the 
near-threshold data on neutral pion production it was recognized that the ppn^ reaction near 
threshold is sensitive to short-range mechanisms in the NN-system. Since the main pion 
exchange term, which dominates the charged pion producing channels, is isospin-forbidden, 
Lee and Riska proposed to counter this shortfall by the additional consideration of pair 
diagrams with an exchanged heavy meson (cr, cj) fT^. Likewise the role of pion-rescattering 
only then was treated with more concern, but also in a somewhat controversial manner, since 
field theoretical models and chiral perturbation theories found different relative signs in the 
TT-exchange amplitude with respect to the direct (Born) term ^3]; 13 • Attempts have 
been pursued to additionally include A(1232) isobar and Su and D13 nucleon resonance 
excitations [ 15-17 ]. To date, calculations in the framework of the Jiilich meson exchange 
model J7], which incorporate all the basic terms like heavy meson exchange (HME), 
delta resonance and off shell effects as well as realistic final-state interactions, yield the best 
results; with the exception of the strength parameter of the HME term the model needs no 
adjustable parameters. 

Close to threshold only a few partial waves have to be considered with Ss, Sp, Ps and Pp 
being the leading ones, whereas the role of / = 2 contributions (Sd and Ds) is still under 
debate. Here the Rosenfeld notation Lplg has been used [12] with Lp being the orbital 
angular momentum of the NN-pair, Iq that of the pion with respect to this pair. In the ppn^ 
reaction the partial wave Sp which apart from Ss governs the isoscalar reaction channels 
dn'^ and pnvr"*", is forbidden by conservation laws. As a result purely isotropic angular 
distributions are found for the ppir^ reaction in the region up to 10 MeV above threshold 
Also, sizeable analyzing powers are only observed for this reaction for excess energies 
Q higher than 16 MeV ^ which follows from the fact that interference effects are only 
possible between amplitudes (Ps, Pp) on the one hand or (Ss,Sd,Ds) on the other. Except 
for Ss, however, the magnitude of all these amplitudes is negligibly small up to Q =16 MeV. 
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For comparison, angular distributions as well as analyzing powers extracted from a study 
of the vr^-producing reactions show clear evidence for the presence of Sp, and hence / = 1 
contributions at energies as low as 2 MeV above threshold [ 10-11 ]. 

Measurements very close to threshold, where all escaping particles are confined to a narrow 
cone around the beam axis, require a perfectly matched detector system which has to cover 
as much of the available phase space as possible in order to keep necessary acceptance 
corrections to a minimum. In Jiilich, the COSY-TOF spectrometer, a scintillator hodoscope 
of cylindrical shape with full 0-symmetry fulfills these requirements to a very large extent. 
Charged particles are detected at polar angles as close as one degree to the beam. In 
the past we have used this spectrometer in meson- and hyperon-production studies and 
pp bremsstrahlung experiments. Parallel to the analysis of the bremsstrahlung reaction 
j2Uj . data obtained simultaneously for pp — > pprc^ were analysed as a cross check. We 
systematically found larger cross sections than given in [T] , j3] . Since the results we deduced 
from the data for the concurring charged vr-producing reactions agreed nicely with the 
published values jTj, jTU], [U]; whereas the cross sections for vr^-production remained high, 
we decided to perform a second experiment to clarify that matter. 

In this paper we will present data from a kinematically complete measurement of the re- 
action pp — > ppiT^ for three energies in the range 292 - 298 MeV, corresponding to beam 
momenta of 796-805 MeV/c. Total cross sections as well as angular and energy distribu- 
tions will be shown and compared to Monte Carlo simulations which are based on phase 
space distributed events that are modified by the final-state interaction (FSI). Results of 
the concurring reactions pp — » c/vr^ and pp pnvr^ which were obtained in parallel will 
only be given to underline the performance of the spectrometer and to prove the consis- 
tency of the analysis. Apart from the excess energy Q we will also use the dimensionless 
parameter t] for labelling energies above threshold defined as rj = qmax/fn-K with q being 
the CM-momentum of the pion and its mass. With p we denote the CM-momentum 
of one proton in the two-proton subsystem. The angles Qq and Bp (in the CM-system) are 
taken with respect to the beam. The momenta q and p are linked by energy conservation 
via y/s = ^/ml + q^ -|- M|p -|- (f- and M^^ = 4(mp -|- p^) where s is the square of the total 
CM energy and Mpp the invariant mass of the two protons. 
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II. Experimental Procedure 
A. Apparatus 

The experiment was carried out with the time-of-flight spectrometer COSY-TOF [21] set 
up at an external beamhne of the 2.5 GeV proton synchrotron COSY at Jiihch/Germany. 
The goal was to study simultaneously and with the same spectrometer all five reaction 
channels which are open in pp-collisions at beam energies near 300 MeV. Besides pp-elastic 
scattering, which is by far the strongest, and the very weak pp bremsstrahlung, our interest 
mainly aimed at the vr-producing reactions pp dn^, pp — > ppn^, and pp — > pnn'^, each 
one with cross sections rising fast with excess energy. 



A schematic sketch of the spectrometer showing the locations of the main components, the 
barrel (B), ring (R), and quirl (Q) detector, as well as the target region is given in fig. 1. To 
prevent the reaction products from undergoing secondary reactions on their way from the 
target to the different detector modules of the spectrometer it is housed within a 3.3 m x 
4 m steel tank which can be evacuated to a pressure of 0.2 Pa. Protons were extracted out 
of the COSY ring in a slow extraction mode allowing spill lengths of up to 10 min. Their 
intensity typically was of the order of several 10^/s. The beam was focussed onto the liquid 
hydrogen target with dimensions 6 mm 0x4 mm |23j. Its front and rear ends were closed 
by 0.9 fim thin hostaphan foils. A set of scintillator veto detectors with central holes of 
various sizes were located 260, 51 and 3.5 cm upstream of the target. They helped to define 
the beam spot in the center of the target to an area with diameter d=3.0 mm. Charged 
particles emitted from the target into the forward hemisphere had to transverse a 0.5 mm 
thin plastic scintillator device f24j which served as start detector. This detector shown 
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Figure 1: Sketch of the COSY-TOF spectrometer. 
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in fig. 2 consists of two concentric rings made out of 16 trapezoidally shaped scintillators 
each, which were placed 30 mm (Ring B) and 50 mm (Ring A) behind the target. The 
inner detector had a central hole of diameter dj=3.3 mm and extended out for 26 mm, the 
outer one had dimensions dj=16 mm and da=130 mm. The minimum angle for observable 
ejectiles was 1.9° if they originated from the center of the target, but smaller values far less 
than 1.0° were possible if the initial pp-reaction took place off the beam axis. The elements 
of the outer ring were arranged in such a way as to show a tiny overlap along their sides 
thus guaranteeing full 0-coverage whereas the ones of the inner ring were only allowed to 
touch each other thereby leaving narrow gaps of order 50 /xm which caused an overall loss 
in 0-coverage of 3 % distributed evenly over the full circle. 




Figure 2: Front and side view of the start detector. The target is placed to the left of ring 
B. 

After a flight path of up to 3 m in vacuum, the reaction products hit the 3-component (B, 
R, and Q) scintillator hodoscope. The circular, 3-layer quirl and ring detector components 
of the endcap were both built according to the same layout. 5 mm thin scintillator sheets 
(BC408) were used as basic material. The first layer in the quirl is made from 48 wedge-like 
scintillators, followed by two layers with 24 left and 24 right wound elements, each element 
cut along an Archimedian spiral In central projection a net-like pattern is formed 

consisting of more than 1000 triangular pixels. A central hole of radius r=4.2 cm was left 
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for the beam, the outer radius of the quirl is 58 cm. The three layers in the ring were made 
from 96 wedges and 2 x 48 left and right wound Archimedian spirals, the inner (outer) radius 
is 56.8 (154) cm, respectively. The barrel consists of 96 scintillator bars, mounted to the 
inside of the tank. Each bar of 2.85 m length has a cross section with measures thickness 
X width =15 mm x 96 mm and is read out from both ends |22j . 

In the course of the experiment which went on for almost two years two slightly different 
set-ups were employed. The one described above with barrel, ring and quirl was used at 
T=293.5 MeV. Azimuthal coverage was complete, the acceptance in polar angle 9 ranged 
from 1.5° — 78°. In an earlier experiment performed at T=292.2 and 298.1 MeV the ring 
component was still under construction, leaving a gap in polar angle 6 between 10.4° — 26°. 
For this run a modular neutron detector [23] was additionally set up downstream around 
the beam line and outside of the tank. It covered approximately the same area as the quirl 
detector and was employed in the study of the concurring pp — * pmT~^ reaction. As a first 
level trigger usually one, in some cases two charged hits were required in the start detector, 
at least two in any of the stop detectors B, R, Q, with 2B, 2R, 2Q and 2(BRQ) denoting 
separately handled trigger patterns. For each of these patterns slightly different deadtimes 
were observed and the necessary corrections had to be applied accordingly. 



B. Performance Checks and Luminosity Monitor 

Simultaneously with the vr-producing reactions, elastic scattering events, which were de- 
tected either as barrel-barrel (BB) or as barrel-ring (BR) coincidences, were written on 
tape. In addition to serving as a tool for calibrating the whole detector system, this re- 
action was used as a luminosity monitor through comparison of experimentally deduced 
elastic scattering events with data taken from the SAID database PHI (see fig. 3). The cross 
section around T = 300 MeV is known with an error of less than 5 %. Any cross section 
of interest determined relative to the one for ppeiast will thus be measured with comparable 
accuracy, since all uncertainties resulting from instabilities in beam intensity and target 
thickness cancel. 

We obtained complementary information on the detector performance from a careful study 
of the pp —>■ (iTT^-reaction. Its cross section at 295 MeV is close to 45 fih fU]. While 
the deuterons with a maximum polar angle at 3° will only reach the quirl, the pions with 
Qmax _ can be observed in all of the three stop-detector components. The unique 
kinematics of this 2-particle reaction allowed the extraction of rather clean and almost 
background free spectra as e.g., invariant d- and vr'^-mass distributions. In the present 
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Figure 3: Measured angular distribution of elastic pp scattering events (solid dots) compared 
to the SAID solution (full line, see text). 



case this reaction was used for an accurate determination of the beam energy Tp, which is 
important in view of the strong Tp dependence of the cross section as well as available phase 
space coverage, which in turn has a strong impact on the acceptance correction. In fig. 4 the 
transverse 7r-momentum pt is plotted versus the longitudinal vr-momentum pi. The boxes 
represent the experimental data, their sizes being proportional to the respective number of 
counts. Concentration of yield is observed along an elliptical curve which aligns nicely with 
the one resulting from two-body kinematics for a beam energy of Tp =293.5 MeV (dashed 
curve). The two solid lines calculated for beam energies which are 0.6 MeV higher and lower 
clearly miss the main body of the data from which we deduce an error in beam energy of 
order 0.3 MeV. 

The modulation in yield along the elliptical curve and the absence of yield for values below 
17 MeV/c is well understood and can be explained as stemming from particular technical 
details in detector construction as well as the rather low pion momenta: pions emitted at 
angles larger than 26° will be observed in the one-layer barrel section; these arc the events 
with pi up to 52 MeV/c. Those with very low momenta, however, will not be able to reach 
it. The ring section covers the range 10.4° — 26°; here, however, a pion has to reach at least 
the second layer in order to fix its point of impact, hence a decrease in yield is observed for 
pion angles approaching the outer range of the ring. Very forward pions with low pt- values 
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Figure 4: Transverse vs. longitudinal TT-momentum in comparison with expectations from 
two-body kinematics for three different beam energies (see text for details). 

were not observed, since the coincident deuterons escaped through the central hole in the 
start detector. 

III. Data Analysis and Calibration 

A. Principle of Measurement 

The photomultiplier signals of all detector elements were used for energy-loss as well as 
time-of-flight measurements. To this end, the analog signals were digitized in Fastbus-QDC 
modules, the timing signals were derived in leading-edge type discriminators and used in 
Fastbus TDCs. Both measurements yield the particle velocity (3. From the particle's point 
of impact the flight direction was obtained from which two additional observables, 6 and (p, 
can be derived. In the endcap with its 3-layer quirl and ring components the flight direction 
is deduced by combinatorical means [21] • For the one-layer barrel elements with two-sided 
readout the point of impact is found from the difference in TDC- values of both ends [22] , 
whereas the TDC-sum is a measure of the flight time. The spatial resolution as given from 
the pixel size in the endcap is A6 = 0.5°, A0 = 2°; roughly the same values hold for the 
barrel elements. Typical values obtained for the time resolution of the start-stop system 
were of the order of 0.5 ns (FWHM) which corresponds to about 2.5 % of the average flight 
time of a /3 = 0.5 particle. 

The momentum 4-vector of each detected particle can then be deduced from the measured 
observables /3, 9, and by applying an additional mass hypothesis. This of course has to 
be checked in comparison with results obtained from Monte Carlo simulations. In case of 
a reaction with one neutral particle like pp ppn^ the 4-vector of that particle is deduced 
by employing energy and momentum conservation and calculating its mass. 
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B. Calibration 

In order to determine 9, (j) and (3 for each particle, several steps in calibrating the detector 
had to be performed After a simple pedestal subtraction in the QDC spectra, a walk 
correction had to be applied for each TDC entry. Depending on the pulse height of the 
scintillator signal, time shifts of up to 3 ns were observed. To correct for these shifts a 
five-parameter fit function was found which relates the observed walk with the pulse height 
given by the corrected QDC value. The complete parameter list for up to 600 scintillator 
channels was determined before each experiment. This was routinely done by use of a laser- 
based calibration system where UV-light from a A^2-laser is fed via quartz fibers into each 
of the individual scintillator channels. Through use of a set of filters the amount of light 
going into the various elements can be varied such as to cover the complete dynamic range 
of the multipliers. 

A check on the differential nonlinearity in the TDC modules revealed no deviations from the 
precision given by the manufacturer (< 0.1 % of full scale), in some cases, however, a channel 
width of 90 ps/channel was observed instead of the preset value of 100 ps/channel. Due to 
different cable lengths and varying transit times in the photomultiplier tubes and electronic 
modules, variations in the signal arrival time arise which are specified as TDC-offsets. These 
offsets were determined in an iterative manner by first comparing overlapping, but otherwise 
identical channels as e.g., neighboring elements in the outer start detector (sect. 11. A) or 
overlap regions from left- and right-wound elements in the endcap. Two-particle reactions 
like PPeiast and dTT~^ with their unique kinematics then were used for the adjustment of all 
remaining detector components. As a byproduct the light velocities within the differently 
shaped scintillator elements were deduced. Whereas in bulk material it is given by c/n (n: 
index of refraction), smaller values were found for the rather fiat, but long elements due to 
an effective lengthening of the lightpath. For the straight quirl elements of wedgelike shape 
and a thickness of 5 mm we found 18.0 cm/ns, for the corresponding ones of the ring 15.8 
cm/ns. A value of 16.1 cm/ns was observed for the 3 m barrel elements. The calibration of 
the neutron detector with its 10 cm thick bars was performed as outlined in 123- Special 
care was applied in the determination of the energy-dependent efficiency which, for neutron 
energies larger than 60 MeV, was found to be rather constant around 13.5 %. 
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C. Monte Carlo simulation 

Throughout the analysis each step was compared to the results deduced from the Monte 
Carlo simulation. The program package written in C++ was developed in Bochum [2H] , [211 
with the purpose of simulating the spectrometer response, thereby giving insight into the 
overall performance, and to determine the detector acceptance. At its core is the CERNLIB 
random event generator GENBOD 32 , which generates a preset number of N-body events 
for a given reaction specified by N, the number of ejectiles, their masses and the 4-momentum 
of the incoming beam. It returns momentum 4-vectors for each particle in the overall center 
of mass system, and weight factors We based on the phase space density of the reaction. 
These weight factors have to be modified, if angular momentum or final-state interaction 
effects are believed to play a substantial role. The Monte Carlo code then boosts this event 
into the laboratory system and each particle is tracked through the complete setup. Nuclear 
interactions are treated according to the INC code [201 which has its roots in the Bertini 
intranuclear cascade model of the late sixties [3Tj. This model assumes the nucleons within 
their respective nuclei to behave like a fermigas and the interactions between the incident 
particle and the nucleus is simulated as the output of an intranuclear particle shower with 
nucleon-nucleon and pion-nucleon reactions. Elastic and inelastic NN and vrN cross sections 
are supplied via extensive data tables. Great care was put into modelling all of the detector 
components in their particular shape as close as possible to reality, including the left- and 
right-wound spirals in the endcap and the trapezoidal elements of the slightly conical barrel 
detector. Options were built into the program to study in detail the influence of the extended 
target and the transverse beam spread or the variation in acceptance due to the narrow gaps 
between neighboring scintillator elements or the size of the central hole in the inner start 
detector. 

In the course of the simulation, interactions of the particles on their way through the various 
active or passive detector parts can be switched on or off, their energy loss and small- 
angle scattering can be investigated separately. Effective low-energy thresholds can thus be 
studied in detail as eg., in case of the endcap, where a hit in at least two layers is needed for 
a successful pixel reconstruction. Also included in the Monte Carlo package is the possibility 
of particle decay. In the present case it was found that a sizeable fraction of the vr """-mesons 
with low momentum had decayed before the detector was reached and a decay myon was 
detected instead. 
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IV. Results 

A. Total cross sections 

By use of energy and momentum conservation, the pp7r° events were identified by calculating 
the missing mass m-^ of the unobserved neutral pion according to 

= (P, + P,-Pi-P2)2, 

where Pb(Pt) denotes the beam (target) momentum 4-vector, respectively, and Pj the one 
of the i-th proton. The high energy part of the distribution obtained at T= 293.5 MeV is 
given in fig. 5 showing a clear signal at the expected location nix — "^n— 0.135 GeV/c^. 
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Figure 5: Missing-mass spectrum of reconstructed pp-K^ -events (dots with error bars) in 
comparison with the Monte Carlo simulation (histogram). The dashed histogram shows the 
normalized spectrum obtained with an empty target. 

The only cuts applied in the analysis were upper limits on the /3-values of the protons in 
order to eliminate events where a much faster pion was wrongly interpreted as a proton. 
Background due to reactions on the target foils or condensates thereon is negligible as 
can be seen from the spectrum obtained from an empty target run and shown by the 
dashed histogram in fig. 5. Also displayed is the result of our Monte Carlo simulation 
(solid histogram) which is in excellent agreement with the data. The angular distribution 
of the inelastically scattered protons is shown in fig. 6 (solid dots with error bars) together 
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with results of our Monte Carlo simulations. A shift towards smaller proton angles and 
away from purely phase space distributed events (dashed histogram) is readily seen. When 
incorporating final state interaction effects into the Monte Carlo code, however, and using 
standard values for effective range Tq and scattering length a, the simulated curve nicely 
reproduces the data as shown by the solid histogram. For more details on our treating of 
the final state interaction see sect. IV. B. 



^ 200 




4 8 12 16 

e [deg] 

Figure 6: Proton angular distribution obtained at T— 293.5 MeV in the laboratory system. 
The solid (dashed) histogram shows Monte Carlo generated results, where final state inter- 
action effects have been included (omitted). 

In a next step total cross sections for the reaction pp ppir^ were extracted from the data 
by summing up the events for m-^ > 0.13GeV/c^ (see fig. 5). By making minute corrections 
for falsely interpreted pp — > pn7r"'"-events of order 3 % from Monte Carlo considerations, 

one finds the number of ppn^-events as given in column 3 of table 1. Column 4 lists the 
acceptances as derived from our Monte Carlo simulations. The main deviation from complete 
acceptance {ace = 1) results from geometrical properties as e.g., the central hole in the inner 
start detector (mainly for the low-energy measurements) and the missing ring component 
(see sect. II. A) for the measurement at T=298.1 MeV. Nevertheless, in all cases acceptance 
values of order 50 % or higher were found. The error in acc was estimated by varying in the 
Monte Carlo code such values as the diameter of the central hole of the start detector (by 
some 10 %) or the beam energy (by 0.2 MeV). An overall error of Aacc — 3.5% was found. 
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Table 1: Total pp —>■ pprc^ cross sections derived at three beam energies. The column labelled 
acc lists the detector acceptance as determined from our Monte Carlo simulations based on 
FSI modulated phase space distributed events. 



T{MeV) 



V 



n 



acc 



jLdt inb-^) 



292.2 ± 0.3 0.29 1586 0.55 



1.25 



2.31 ± 0.06 ± 0.23 



293.5 ± 0.3 0.30 2524 0.64 



1.41 



2.80 ± 0.06 ± 0.21 



298.1 ± 0.2 0.35 4426 0.49 



2.43 



3.72 ± 0.06 ± 0.28 



The integrated luminosities as given in column 5 of table 1 were derived from a comparison 
with the results obtained for pp-elastic scattering in the angular range 32° — 55°, a range 
that was covered by the barrel alone. For the measurement at T=293.5 MeV, where the 
ring was available as well, the quoted value not only is corroborated by the results found 
for the extended range 18° — 70°, but also from a comparison with results for pp —>■ dn^. 
For this reaction a total cross section of cr = 42 ± 5 /xb was found in very good agreement 
with an interpolated value of 44 fih deduced from For the pp —>■ pnir'^ reaction at 

T=298.1 MeV we determined a cross section of a +=3.5±0.4 fib which agrees nicely with 



the interpolated value of 3.1±0.4/i6 extracted from [7], [S]. 

Total cross sections for the reaction pp ppn^ are given in the last column of table 1 
together with their statistical and systematical errors. The latter are composed of 5 % from 
the pp cross section, 3.5 % from the acceptance correction, 1.5 % from the uncertainty in 
misinterpreted background events and 2 % in the luminosity determination. 
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B. Dalitz Plots and Energy Distribution 

As mentioned before, a crucial point in the determination of the total cross section is the 
extrapolation of the measured data into the full region allowed by phase space. A thorough 
understanding of the detector response is essential, as is the knowledge of any deviation in 
the original (physical) data from those given by simple phase space, i.e. s-wave distributed 
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Figure 7: top: Dalitz plots of experimental (left) and simulated (right) data with kinematical 
limits given by the solid lines. The enhancement due to FSI is seen in the upper left corner of 
the event distribution, bottom: projection of the Dalitz plots onto the M^p-axis. Experimental 
data are given by solid dots with statistical error bars, the solid histogram shows the simulated 
distribution including FSI effects. The dashed histogram represents the simulated distribution 
for pure phase space. 



events. A convenient way to compare in a comprehensive manner experimental and sim- 
ulated data from a 3-particle reaction is to set up a Dalitz plot, which quite instructively 
allows the extraction of the underlying physics from the event distribution. Purely phase 
space distributed events will result in a uniformly covered Dalitz plot, deviations from such 
a distribution will induce some modulations. The Dalitz plot as found from the present 
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(kinematically fitted) data at T=293.5 MeV is shown in fig. 7 and compared to the result 
of our Monte Carlo simulation (top row left and right, respectively). As a consequence of 
the strong final state interaction in the pp system a systematic enhancement in the top left 
corner is observed in the data. The simulated distribution was obtained by incorporating 
the FSI formalism as developed in [HSI and later refined by Morton et al. 36j. In short, 
one calculates additional weight factors Wfsi which are multiplied with the event weight We 
given by the CERNLIB subroutine GENBOD (see also sect. III.C). These factors are given 
as 

= . [C- . TCM + Mi ^^r, . TCM _ ^ _ ■ ^7.)) V^ 

where T^J^ denotes the pp center of mass kinetic energy = Mpp — 2mp and the 

Coulomb penetration factor 

— e^^-fP-l 

with 7p = ""'f^pp'^ , Here a is the fine structure constant, Ppp = ^j2fippT^p^ and fipp is the 
reduced mass of the pp-system. The term with /i(7p) = 7p ■ J2n n(n'^+-y^) ~ 0-5772 — Inyp 
contributes only little and is often omitted. From literature we took the standard values 
ao=-7.83 fm and ro=2.8 fm [37j as input parameters for the scattering length and effective 
range, respectively, for the two protons in the ^5*0 state. The smooth rise in yield which 
readily is seen in the simulated data is partially obscured in the experimental data due to the 
limited statistical accuracy. However, when inspecting the projection of the 2-dimensional 
Dalitz plot distribution onto the x-axis it is obvious that only through inclusion of FSI effects 
the simulated data follow the experimental ones (fig. 7, bottom) whereas a simulation based 
on purely phase space distributed events fails over the whole range. 

The proton energy distribution T^^^ at 293.5 MeV in the center of mass system is shown in 
fig. 8. That close to threshold proton energies only range up to roughly 7 MeV with a very 
distinct maximum near 1 MeV. The FSI-based simulation given by the dashed histogram 
reproduces the data very well, whereas a simple phase space distribution (given by the 
dotted histogram and normalized to the experimental yield) fails completely. Also shown is 
the result of a simulation where final state interaction effects have been included but with 
a much reduced scattering length of ao=-1.5 fm as was suggested in ^ (fig. 8, dash-dotted 
histogram). When normalizing this distribution to our experimental data as was done in 
fig. 8 very poor agreement is observed. If on the other hand one normalizes only to the 
upper part beyond T^^^=3.0 MeV the maximum around 1.0 MeV is totally missed and 



15 



0123456789 10 

[ MeV] 

Figure 8: Proton energy distribution T^J^^ in the 2-proton CM system (solid dots with statis- 
tical error bars ) in comparison with the Monte Carlo simulation including FSI with standard 
parameters (dashed histogram). The extrapolation to full acceptance is shown by the solid 
histogram. The pure phase space distribution is given by the dotted histogram, a simulation 
with FSI and a scattering length ao = — 1.5/m by the dash-dotted histogram. The insert 
shows Qp"^^ , the smaller of the two proton angles, vs. T^J^^ (see text). The solid line denotes 
the upper boundary in Q^^^ as deduced from our MC simulation, the lower one coincides 
with the X-axis showing clearly the limit in our acceptance slightly below 2°. 

we cannot see how a proper extrapolation can be achieved. From the insert showing in a 
two-dimensional plot Qp^^, the smaller of the two proton angles, versus T^^^, it can be 
seen, that the maximum yield found near T^^^=l MeV is connected with one of the two 
protons emerging towards 0- values well below 5°, a range that is almost completely without 
detector coverage in each of the two internal target experiments. In view of this observation 
we argue that a description of the data and hence an extrapolation towards full acceptance 
can only be performed correctly by using the standard values ao=-7.83 fm and ro=2.8 fm as 
scattering length and effective range, respectively. The result is given by the solid histogram 
in fig. 8. 

One important remark should be added. In the work of Flammang et al. on the pp —>■ pmr'^ 
reaction it was argued that the acceptance calculation and thus the extraction of a total 
cross section is insensitive to the presence, absence or details of the FSI model ^8J, since 
the additional weight factors used for the extrapolation cancel out. This statement is true 
to some extent and has been checked in the present analysis; the extrapolation to full 
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acceptance as given in fig. 8 by the solid line could also be obtained from our experimental 
points by neglecting in the Monte Carlo simulation any FSI effects and assuming purely 
phase space distributed events. The agreement between experimental and simulated data 
in that case, however, would have been of that quality shown in fig. 8 by the dotted line. 
Instead, and we like to stress this as a crucial point, our initial distribution, although 
modified by the detector acceptance, not only showed the enhancement due to FSI, but this 
distribution was also well reproduced by our Monte Carlo simulation, which we had based 
on i) our knowledge of detector performance and acceptance and ii) the use of standard 
FSI parameters for scattering length and effective range. If we would have started from 
a distribution with (almost) no enhancement nothing would have been gained from an 
extrapolation where one tries to adjust the Wfsi weights alone. 

C. Angular Distribution 

Acceptance corrected angular distributions for the pion- and relative two-proton-momentum 
vectors are shown in fig. 9 for the measurement at T=293.5 MeV. The solid lines are the 
result of fits in terms of Legendre polynomials Pi to the experimental data W{cosQ) oc 
1 + ' Pi{cosQ) up to Z = 2. Only very small a2-coefficients were found as expected 

that close to threshold, namely 02 = +0.05 ± 0.10 for the top and 02 = +0.08 ± 0.11 for the 
bottom distribution. Similar results were obtained at the other two bombarding energies. 
The observed symmetry around cos G = (or the vanishing of the Pi-term) which stems 
from the entrance channel consisting of two identical particles was not used as a constraint in 
the fit; instead we consider it as a valuable check on the quality of our acceptance corrections. 
The angular distribution obtained for the pp —>■ pm^^ reaction at T=298.1 MeV correspond- 
ing to ?7 = 0.19 was found to be consistent with isotropy which is in agreement with the 
published data f5|. Sizeable anisotropics were found for the angular distribution of the 2- 
particle reaction pp —>■ rfvr"^ yielding an a2-coefficient of order +0.2 and larger with, however, 
an error of comparable size. This is a consequence of the significant acceptance corrections 
which had to be performed to account for those pions which had decayed before reaching 
the detector, as well as those which were too low in momentum as being able to reach the 
second layer of the endcap (see also sect. III.C). 



17 



-1 -0.8 -0.6 -0.4 -0.2 0.2 0.4 0.6 0.8 1 

COS 0q 

,— , 0.5 F 

(/I - 

a 0.3 u 

0.1 r 

F 

-1 -0.8 -0.6 -0.4 -0.2 0.2 0.4 0.6 0.8 1 

COS 0p 

Figure 9: Acceptance corrected angular distributions in cos Qq (direction of pion momentum) 
and cosQp (direction of relative two-proton momentum) in the CM system as obtained at 
T=293.5 MeV. The solid lines are the result of Legendre polynomial fits up to I = 2. 

V. Discussion 

The total cross sections listed in table 1 are compared to results found in the literature 
[ 1-3 ], [ 5-6 ] in fig. 10. In the range rj = 0.2 — 1.0 a rises over two orders of magnitude up 
to values close to 100 /ib. The present data exceed the ones from lUCF P and CELSIUS [S] 
by roughly 50 %. We connect this discrepancy with the one basic difference in experimental 
set-up, namely theirs being an internal, ours, however, an external target experiment. The 
most critical step in extracting an absolute cross section lies in the acceptance correction, 
which depends on the knowledge of i) geometrical coverage of the available phase space, 
ii) angular distribution effects, and iii) deviations from s-wave distributed events due to 
final state effects. Since, most probably, the geometrical coverage is well accounted for in 
all three cases, and angular distribution effects are negligible, we believe an underestimate 
of the strong FSI effects to be responsible for the large discrepancy. Our reasoning is as 
follows: The lower the energy above threshold, the more important the relative weight of 
FSI becomes, and, simultaneously, the more yield gets compressed into the very small angle 
region. If on the other hand the minimum angle is only near 4°, as is the case in the 
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internal target experiments at lUCF P] or CELSIUS [3], more and more yield is missed. 
To extrapolate into the uncovered region requires a very profound knowledge of the original 
distribution, and most often some model dependence will enter the extrapolation procedure. 
As was mentioned before the data given in jl] were corrected by using non-standard FSI 
parameters (scattering length ao=-1.5 fm), an approach which completely failed to reproduce 
our data. The low-energy CELSIUS-data [H] were obtained by detecting the two decay 
photons instead of the two emerging protons, a method with an acceptance which is only 
weakly dependent on excess energy and is free of any threshold effects, but suffers from the 
very small geometric coverage of only 5 %. A rather elaborate acceptance correction has 
been developed for the high-energy CELSIUS data [H] where FSI- and angular distribution 
effects were considered. As can be seen from fig. 10 these data exceed the previous lUCF 
data jjy as well although by a smaller margin of some 15%. 
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Figure 10: Summary of cross section data as known today forpp ppir^ up to rj = 1.0. The 
present data given by solid dots clearly are higher than the earlier ones from lUCF (open 
circles, Jl^) and CELSIUS (solid downward triangles, ^3]). 

A summary of Legendre polynomial expansion coefficients 02 found in the literature for the 
pion producing reactions in the near threshold region is shown in fig. 11 as a function of f]. 
Three bands of data are discernible. The fastest rise with rj is observed for the pp diT~^ 
reaction (given by open symbols, refs. [10 — 11], ^38 ) where deviations from 02 = are 
already found below rj = 0.1. The positive a2-coefficients are an indication for strong p- 
wave contributions. Much less information is available for the other vr^ -producing reaction 
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Figure 11: Comparison of experimental a2- coefficients as obtained from Leg endre polynomial 
fits to 7T-angular distributions in the near threshold region, pp dn'^-data are given by open 
symbols, the ones for pp —>■ pmc'^ by solid asterisks and diamonds. Solid dots, squares and 
triangles as well as open asterisks are used to denote the data forpp ppir^ (see text). 



measured at lUCF and TRIUMF and plotted as asterisks [S] and diamonds , respectively. 
Yet a comparable behaviour of a fast rise can be observed as well, which, however, only starts 
somewhat later near r] = 0.2. For both these reactions with an isoscalar two-nucleon pair 
in the exit channel the Sp partial wave which is not Pauli forbidden contributes almost 
from threshold. A totally different pattern is observed in case of the pp ppvr "-react ion. 
Near-isotropy is found for almost the whole //-range shown in fig. 11 with 02 staggering 
between -0.2 and -1-0.2. It is not clear whether d-wave contributions which induce negative 
02 coefficients might be present that close to threshold. A unique trend (toward positive 
values) is only observed for rj > 0.9. 

To get a more quantitative estimate on the amount of p-wave contributions to the ppn^- 
reaction at 77 = 0.3 we used the fact that close to threshold the dynamics for pion production 
as given by the matrix element is almost energy independent. Instead the rapid rise in 
cross section is governed by kinematical effects since phase space as well as the radial wave 
functions of the participating particles change drastically. In a non-relativistic treatment the 
available phase space volume scales as q-p. The wave functions can be approximated by their 
asymptotic form which is valid for small arguments, and from this one obtains a factor q'"' -p^^ 
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Table 2: Reduced values derived from fitting Monte Carlo generated momentum distribu- 
tions to our experimental data by varying the amount of Ps or Pp contributions (see text). 



Lplq 


Ps 


Pp 


f 


0.05 0.10 0.20 
1.24 1.41 1.93 4.15 


0.05 0.10 0.20 
1.24 1.26 1.36 1.90 



with Iq (Ip) denoting the angular momentum of the pion (the NN pair), respectively. Hence 
each of these additional terms has a particular momentum dependence as e.g., aps (x q ■ p^ 
and app oc -p^ To account for these higher order terms we introduced into our Monte 
Carlo code additional weight factors wps = Cp^ ■ qp^ and wpp = Cpp -p^q^ with adjustable 
fit parameters Cp^ and C|,p, which were added to the dominant weight We ■ Wfsi discussed 
in section IVB. Due to the limited statistical accuracy in our data we did not attempt a 
combined fit where the Ps and Pp fractions were varied simultaneously. Instead we fitted 
our experimental q- and p-momentum distributions by a Monte Carlo generated sum of 
either Ss+Ps or Ss+Pp terms. In view of the good overall acceptance of our detector with 
only a small central hole of slightly less than 2° any acceptance dependence of the weight 
factors, however, was neglected. The resulting x^-values as obtained for various fractions 
f are shown in table 2 where f is given through aps{p) = f ■ a and ass = (1 — /) ■ c"- The 
best fit with a reduced x"^/^ = 1-24 with n=65 degrees of freedom is found for pure s-wave 
behaviour. Up to 10% of a Pp contribution is allowed by the fit, Ps contributions are even 
weaker. As a result we assign to the Ss-strength for the ppvr*^ reaction at ?7=0.3 a fraction 
(l-f)=0.95±0.05, leaving only a minute contribution of 0.05±0.05 for Pp and Ps waves. 

In fig. 12 we show the presently deduced partial cross sections as a function of r] together 
with data found in the literature [ 4-6 ]. Although some discrepancies in the data exist, 
which we attribute to the different methods applied, the general trend seems to be pretty 
much established. Up to ?7=0.7 Ss given by the diagonally hatched band is strongest, but 
gets almost negligible above ?7=0.9. At the same time the Pp strength (cross-hatched band) 
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starts to dominate, but saturates near 0.6. Ps strength (vertically hatched band) seems 
to be present already below 7^=0.4 and exceeds Ss near ?7=0.8. In the crossover region the 
discrepancies are the largest which most probably is due to the different methods used. 
The studies carried out at CELSIUS and COSY-GEM were also performed with 
unpolarized beam and target; from these experiments sizeable Pp strength was reported 
well below ri=0.6. A much higher sensitivity can be expected when measuring polarization 
observables as was recently done at lUCF In that study aps was found to be larger than 
app near 7^=0.56. 




0.2 0.4 0.6 0.8 1 1.2 

r[ 

Figure 12: Summary of partial cross sections ctli/<^ vs. rj. The present results are given 
by the full circles. Solid symbols (squares and triangles) which form the cross-hatched band 
denote app, open squares and triangles in the vertically hatched band represent apg. ass is 
given by stars and asterisks (diagonally hatched band). 

Up to ?7 = 0.7 the Ss partial wave dominates the ppvr^-reaction. As was pointed out in 
|17j the cross section in this region could only be reproduced satisfactorily when adding an 
extra heavy meson exchange (HME) contribution with a "free" strength parameter to the 
coherent sum of direct and rescattering term, both of which make up roughly 60 % of the 
measured yield. Yet at energies above rj = 1 where even higher partial waves with I > 2 
start to contribute significantly and where A-isobar contributions have to be considered the 
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predictions underestimate the data more and more. In view of the present findings of a 
much larger cross section also at low t] values which would require an even larger strength 
parameter for HME contributions it might be conceivable that something more essential is 
missing in the theoretical description of the near-threshold ppir^ reaction. It seems to be 
worth mentioning that poor agreement is also found in recent calculations of polarization 
observables [iU] . 

VI. Summary 

The pp —>■ ppiT^ reaction has been investigated in a kinematically complete experiment by 
detecting the two protons in the large acceptance spectrometer COSY-TOF set up on an 
external beam line of the proton synchrotron COSY. The calibration of the detector and 
the luminosity determination was carried out by analysing the pp-elastic scattering reaction 
recorded simultaneously. Extensive Monte Carlo simulations were carried out in order to 
determine the overall acceptance with high precision. Results derived for the concurring 
reactions pp dn^ and pp —>■ pmr^ are in very good agreement with data from the 
literature. Total cross sections for the reaction pp ppn^ were obtained at bombarding 
energies of 292.2, 293.5 and 298.1 MeV, which exceed previous measurements by roughly 
50 %. Invariant mass and energy distributions are governed by strong final state effects. The 
importance of detecting protons well below 4° to the beam axis is stressed. The angular 
distributions for the pion and relative two-proton momentum vectors were measured for 
r^-values as low as ?7 = 0.3 and were found to be nearly isotropic as expected that close to 
threshold. From an analysis of the corresponding momentum distributions the combined 
fraction for Ps and Pp partial cross sections was estimated to be 5 ± 5 %. 

Acknowledgements 

The help of the COSY crew in delivering a low-emittance proton beam is gratefully ac- 
knowledged. The data are based in part on the analysis work performed by B. Jakob fI7\ . 
We like to thank R. Klein and M. Wiirschig-Porsel for their continued assistance in solving 
technical problems. Helpful discussions with C. Hanhart are very much appreciated. Fi- 
nancial support was granted by the FFE fund of the Forschungszentrum Jiilich and by the 
German BMBF. 



23 



References 

[1] H.O. Meyer et al, Nucl. Phys. A539 (1992) 633 

[2] S. Stanislaus et al., Phys. Rev. C44 (1991) 2287 

[3] A. Bondar et al., Phys. Lett. B356 (1995) 8 

[4] H.O. Meyer et al., Phys. Rev. C63 (2001) 064002 

[5] R. Bilger et al., Nucl. Phys. A693 (2001) 633 

[6] M. Betigeri et al. (GEM collaboration), Phys. Rev. C65 (2002) 064001 

[7] J.G. Hardie et al., Phys. Rev. C56 (1997) 20 

[8] R.W. Flammang et al., Phys. Rev. C58 (1998) 916 

[9] W.W. Daehnick et al., Phys. Rev. C65 (2002) 024003 

[10] M. Drochner et al., Phys. Rev. Lett. 77 (1996) 454 

[11] P. Heimberg et al., Phys. Rev. Lett. 77 (1996) 1012 

[12] T.-S. Lee and D.O. Riska, Phys. Rev. Lett. 70 (1993) 2237 

[13] B.-Y. Park et al., Phys. Rev. C53 (1996) 1519 

[14] T.D. Cohen et al., Phys. Rev. C53 (1996) 2661 

[15] U. Van Kolck et al., Phys. Lett. B388 (1996) 679 

[16] M.T. Pefia et al., Phys. Rev. C60 (1999) 045201 

[17] C. Hanhart et al., Phys. Lett. B444 (1998) 25 

[18] C. Hanhart et al., Phys. Rev. C61 (2000) 064008 

[19] A.H. Rosenfeld, Phys. Rev. 96 (1954) 130 

[20] COSY-TOF collaboration, Phys. Lett. B429 (1998) 195 

[21] M. Dahmen et al., Nucl. Inst, and Methods in Phys. Res. A348 (1994) 97 

24 



A. Bohm et al., Nucl. Inst, and Methods in Phys. Res. A443 (2000) 238 

A. Hassan et al., Nucl. Inst, and Methods in Phys. Res. A425 (1999) 403 

P. Michel et al., Nucl. Inst, and Methods in Phys. Res. A408 (1998) 453 

L. Karsch et al, Nucl. Inst, and Methods in Phys. Res. A460 (2001) 362 

R.A. Arndt et al., Phys. Rev. C62 (2000) 034005; see also: Scattering Analysis Inter- 



active Dial-in, URL: http: //said. phys. vt.edu 



B. Jakob, Ph.D. Thesis, Technische Universitat Dresden (2001) 
S. Brand, Ph.D. Thesis, Ruhr- Universitat Bochum (1995) 
U. Zielinsky, Ph.D. Thesis, Ruhr-Universitat Bochum (1999) 

P. Cloth, in MC93, Int. Conf. on MC Simulations in High Energy and Nuclear Physics, 
edited by P. Dragovitch, S.L. Linn, M. Burbank (World Scientific, Singapore 1994), 
p. 56 

H.W. Bertini, Phys. Rev. 188 (1969) 1711 

GENBOD, CERN Program Library Long Write-up W515 (1993) 
M. Drochner et al., Nucl. Phys. A643 (1998) 55 
K.M. Watson, Phys. Rev. 88 (1952) 1163 

A. B. Migdal, Sov. Phys. JETP 1 (1955) 2 

B. J. Morton et al., Phys. Rev. 169 (1968) 825 

H.P. Noyes, Ann. Rev. Nuclear Science 22 (1972) 465 

B. G. Ritchie et al., Phys. Rev. Letters 66 (1991) 568 
R.G. Pleydon et al., Phys. Rev. C59 (1999) 3208 

C. Hanhart, Proceedings MESON2002, Krakow/Poland (2002), |^^th/0207Q15 



25 



